Atmospheric Dependence of the Stopping Cosmic Ray Muon Rate at Ground Level 
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The rate of low energy (< 150 MeV) cosmic ray muons was measured at ground level as a function 
of several atmospheric parameters. Stopped muons were detected in a plastic scintillator block and 
correlations were determined using a linear regression model. A strong anti-correlation between 
fractional changes in the ground-level pressure and stopping muon rate of —3.0 ± 0.5 was found, 
and also a —4.1 ± 0.5 anti-correlation with the fractional change in atmospheric height at 10 kPa 
pressure. A weak positive correlation with the 10 kPa temperature was also found, but it was shown 
not to be statistically significant in our data set. The same analysis was applied to the total rate 
of all charged cosmic ray particles detected with the same apparatus, and good agreement with 
previous work was seen. The pressure and height correlation parameters for stopping muons are 
larger than for the total rate of all charged particles by factors of about 1.6 and 3.7, respectively. 

PACS numbers: 96.50.S-, 96.50.sb, 94.20wq 
Keywords: cosmic rays, muons, atmosphere 



I. INTRODUCTION 

The rate of cosmic rays at ground level varies in time 
due to changes in both the incident flux at the top of the 
atmosphere and changes in atmospheric conditions. (For 
an overview see Ref. PQ.) Variations in the primary spec- 
trum are thought to be due in part to the indirect effect of 
the magnetospheric distortions caused by solar activity, 
and these "Forbush" variations can occur on a timescale 
of hours and days and be up to about 20% in magni- 
tude. Separately, the evolution of cosmic ray showers 
in the atmosphere is affected by gas density variations, 
which are in turn related to measurable pressures and 
temperatures at various altitudes, as has been discussed 
in the literature for decades. Incident protons between 
roughly 15 and 20 km altitude initiate collisions with nu- 
clei that produce primarily pions and spallation nuclear 
fragments. When this occurs at comparatively higher al- 
titude and at lower air density, the pions are more likely 
to decay to muons rather than to reinteract to create 
lower energy pions and nucleonic secondaries. Higher 
average energy decay muons, though fewer in number, 
are less likely to suffer enough ionization energy loss to 
stop and decay before reaching the ground. On the other 
hand, if the pionic interactions take place at compara- 
tively lower height and higher air density, pions are more 
likely to reinteract before decaying, ultimately produc- 
ing a larger number of lower energy muons closer to the 
ground. The interplay of these processes leads to the 
well-known broad momentum spread of muons and other 
particles at ground level [2J, but also to correlations of 
particle rates and atmospheric conditions. 

The ground-level air pressure and cosmic ray rate are 
anti-correlated, since denser air causes more ionization 
energy loss to occur, and therefore earlier stoppage of 
muons. Also, for a given temperature lapse rate, denser 
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air at ground level means a higher altitude at which 
muons are formed, which in turn gives muons more 
chance to decay before arriving at the surface. 

The mean height of muon formation coincides with 
roughly the 10 kPa (100 mb) pressure level. Again, an 
anti-correlation is expected with changes in this height, 
since formation of muons at higher altitude means that 
fewer will arrive at ground level. Direct measurements of 
the height of the atmosphere using balloon-type sound- 
ings are available j3j 0|. We used the 10 kPa pressure- 
level height to analyze our data. We then studied the 
variability of the results with respect to this choice, as 
will be shown. 

The temperature of the atmosphere at the mean muon 
production altitude may be supposed to affect the muon 
rate in the following way: if the air is warmer it is less 
dense (at fixed pressure) , so the pions have greater prob- 
ability of decaying to muons before they reinteract. This 
would favor creation of higher average energy muons that 
reach ground level, and also lead to larger numbers that 
penetrate deep underground. A positive correlation with 
temperature is seen underground, for example in Refs. 5 
7 . But the expected effect at ground level is less definite, 
since lower average energy muons created from pions in 
a denser atmosphere may be more copious in number 
at the surface, formed together with fewer higher energy 
muons. Indeed the temperature correlation was too small 
to measure in Ref. [5] for vertically incident muons in a 
sea-level telescope of counters. 

Shower creation and propagation can be modeled 19| 
or simulated in detail, including the correlations among 
the various parameters characterizing the atmosphere, 
but this was not the goal of the present study. Our goal 
was the experimental determination of the correlation 
among several atmospheric parameters and the ground- 
level stopping muon rate. We also measured the param- 
eters of the total rate of charged particles to see whether 
they are measurably different. The properties of the at- 
mosphere that we examined in relation to the particle 



rates were as follows. First, we used the ground-level 
air pressure, P, measured at the site of the detector. 
Second, we used weather data giving the height of the 
atmosphere at certain specific pressures, H, via public 
data, [3 . Third, from the same weather data set, we used 
the temperature of the atmosphere, T, at certain specific 
pressures. Finally, we looked for correlations with the 
low-altitude humidity of the air, but these were found to 
be nonexistent in this work. Atmospheric pressure at any 
altitude is proportional to air density, as is the inverse of 
air temperature, and it is the density that ultimately reg- 
ulates the evolution of cosmic ray showers. Our results 
were more satisfactory, however, when looking for the 
separate correlations with the given parameters. 

In order to investigate the rate dependence of cosmic 
rays and stopping muons, the following linear regression 
model was adopted. As discussed in the next section, 
a seven week span of time during which the data did 
not exhibit drastic fluctuations in rate <£, was defined 
by inspection of the data from a five month long run. 
For that span of time the average values of all variables 
in the data were established (bracketed symbols). The 
expression for relating the variables was taken to be 
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where a, j3, and 7 are unitless coefficients of fractional 

change, and the other variables were defined above. This 
expression can be viewed as the linear part of the Taylor 
series that should, to good approximation, capture what- 
ever the effective relationship among the given variables 
is. A least-squares fit was used to determine a, j3, and 
7, and reduced \ 2 g ave a measure of the goodness of the 
fits. The parameters are assumed to be uncorrelated for 
the purpose of comparing our results to historical prece- 
dents for non-stopping muons. Note that in the literature 
there are various names and formulations for these coeffi- 
cients. For example, a is sometimes called the "negative 
partial barometer coefficient" and reformulated in units 
of %/cm Hg. If /? and 7 are not fitted, it is called the 
"total barometer coefficient" . 

Section [Hi will discuss our equipment and procedures 
for these measurements, followed by the results and com- 
parison to other work in Section [III| Conclusions will be 
stated in Section [TV] 



II. APPARATUS AND PROCEDURES 

The detector for both stopping muons and for counting 
the total rate of through-going cosmic rays is shown in 
Fig.l. It consisted of a 30 cm cube of polystyrene scin- 
tillator, viewed from opposite sides by two 5" diameter 
photomultipliers (PMT). For vertically impinging parti- 
cles, this detector stopped muons of kinetic energy less 
than 100 MeV (177 MeV/c), while for oblique tracks the 
maximum energy was 150 MeV (233 MeV/c). This de- 
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FIG. 1. Scintillator block and photomultiplier pair used to 
detect both stopping muons and to count passing cosmic rays. 



fined what we mean by low-energy stopping muons, of 
which we determined the rate variations. There were 
no separate trigger counters to select vertically-moving 
particles, so this apparatus accepted the full-sky angu- 
lar distribution of arriving cosmic rays. The equipment 
was located under a thin metal roof in a building with a 
largely unobstructed view of the full sky. 

The signals from each photomultiplier were sent to 
discriminators which formed logic pulses from all input 
signals larger than a pre-set threshold of 120 mV. The 
discriminator signals were sent to a coincidence unit to 
create pulses only when signals from both tubes were 
present within a 20 nsec time window, thus reducing noise 
and allowing only one pulse per particle. These pulses 
were used for determining the muon lifetimes with use 
of a time to digital converter (TDC). The data were ac- 
cumulated event-by-event using a National Instruments 
data acquisition (DAQ) unit (Model NI USB-6221) and a 
LabVIEW-based control system. A stopping muon can- 
didate was defined by a pulse pair with a TDC "Start" 
pulse followed by any "Stop" pulse within 20 //sec. Pulse 
pairs with a longer time interval were plentiful, caused by 
uncorrelated cosmic ray particles, and these were counted 
but not individually recorded. The PMT singles rates 
were ~ 300 Hz, the coincidence total count rate for the 
experiment averaged 24/sec, while the rate of stopping 
candidates was about 6/min. For each candidate stop- 
ping muon event, the following parameters were recorded: 
a time-stamp, the decay time in //sec with 12.5 nsec (80 
MHz) time resolution, three redundant pressure readings, 
and the total of cosmic ray event pairs detected up to that 
moment. 

The ground-level pressure was measured simultane- 
ously with the cosmic ray data using the same LabVIEW- 
based DAQ system. We used three single-chip barome- 
ters (Motorola MPX4115A) to record the air pressure 
whenever a stopped-muon event occurred. The average 
of the three readings was used as the final value for the 
pressure. We used METAR weather data from the Al- 
legheny County Airport (KAGC) located 16 km away, in 
order to check the accuracy of the lab pressure data. 



The atmospheric data at high altitude was obtained 
using publicly available Integrated Global Radiosonde 
Archive (IGRA) data [4 . A parser of the extensive ASCII 
data set was written to extract the needed pressure, al- 
titude, temperature, and time information. The data 
from the in-lab detectors and the atmospheric data were 
merged together chronologically in the off-line analysis. 
The IGRA data was available in 12 hour intervals. Ini- 
tially we binned the muon data in one-hour intervals to 
get the average pressure and the number of both stop- 
ping and total events. These data were then rebinned, as 
needed, to match the available IGRA atmospheric data. 
Temperature and atmospheric height data were available 
in this data set at only a small set of pressures. One of 
them was 10 kPa, which occurs at about 16.6 km altitude, 
where primary cosmic rays have started to shower. This 
was the pressure level selected by us to fit our model. 

The experiment ran for 21 weeks from July 5, 2011 to 
December 1, 2011, with some gaps for equipment repair. 
A total of 1.31 x 10 6 stopping muons were detected in 
this time, for an average rate of close to 350 per hour. 
At the low muon momenta considered in this measure- 
ment, there is a slight (~ 10%) excess of /i + over fi~ 
particles [2]. Our apparatus did not distinguish between 
the two. Stopped \i~ events are known to undergo atomic 
capture on carbon atoms in the scintillator, leading to a 
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reduced lifetime due to the weak /i 
The value is 2.028 ± 0.002 ^sec [TO]. The fi + events de- 
cay with their free-space lifetime of 2.1969811±0.0000022 
/usee [2J. The measured muon lifetime from the effective 
exponential decay distribution was 2.117 ± 0.003 ^sec, 
between the two other values, as expected. This verified 
that our data sample consisted of muons that stopped 
and decayed. Within the 20 /xsec time window for se- 
lecting stopped muons, the fraction of uncorrelated non- 
stopping background track pairs was 11.5%. This back- 
ground was subtracted in the following stopped muon 
analysis. 



III. RESULTS 

Figure [2] shows the resultant fractional rate variation 
of stopping muon candidates over the entire data period, 
together with the model fit discussed below. The rate is 
presented as the fractional deviation from the mean of 
the fitted period between July 19 to September 8, 2011. 
The average variation during this time is zero, by con- 
struction. The error bars are purely statistical. One sees 
day-to-day fluctuations of about ±2.5%, as well as larger 
excursions outside the fitted region of up to 20%. 

In seeking correlations with atmospheric data, it was 
found that the large jumps visible in the second week 
and again after the eleventh week were impossible to 
match. We ascribed these large jumps to changes in 
"space weather" of the Forbush variety that were out- 
side of our scope to track and model. An interval of data 
that was comparatively devoid of drastic fluctuations in 



TABLE I. Model parameters for the correlation of particle 
rates with ground-level pressure (a), the 10 kPa height of 
the atmosphere (/}), and the 10 kPa temperature (7). The 
third column is for low energy muons stopping at ground level, 
while the fourth column is for the total detected cosmic par- 
ticle rate. 



Parameter 


Eq. 


tl 


Stopped 
Muons 


Total Particles 


Pressure 


Q 




-3.2 ±0.5 


-1.94 ±0.10 


Altitude 


P 




-2.7 ±0.9 


-0.8 ±0.2 


Temperature 


7 




±0.35 ±0.17 


±0.08 ±0.04 








xl = 1-07 


x2= 1-09 



TABLE II. Model parameters for the correlation of particle 
rates as in Table IT] but with only two free parameters. The 
third column is for muons stopping at ground level, while the 
fourth column is for the total detected cosmic particle rate. 



Parameter Eq. 1 Stopped 
Muons 


Total Particles 


Pressure a —3.0 ± 0.5 
Altitude -4.1 ±0.5 


-1.9±0.1 
-1.1 ±0.1 


xl = 1.11 


Xl- 1.13 



rate was used to fit to the model. This was the seven 
week period marked in the figure as "Model Fit Range" . 
The data and the model fit (Eq. fTJ) from this restricted 
range of data are shown in Fig. [3T The solid red curve 
in both figures is the result of the regression model fit. 
Each of the three parameters in the model is a unitless 
scale factor relating a fractional change in atmospheric 
condition to a fractional change in particle rate, as seen 
in Eq. [TJ A least-squares fit optimized parameters a for 
the ground-level pressure, f3 for the 10 kPa atmospheric 
height, and 7 for the 10 kPa temperature. The resul- 
tant values of those parameters are given in Table [TJ The 
uncertainties on the parameters come from the diagonal 
elements of the error matrix associated with the least- 
squares fit. Our data fits the regression model ansatz 
with good agreement, as demonstrated by the reduced 
chi-squared xl value of 1.07, allowing an rms spread of 
y/2/v = 0.15. The Pearson's correlation coefficient is 
0.73, which is significant in view of the large statistical 
uncertainties on the data points. 

The visual appearance of the fit within the fitted range 
is good. We note that there is an anti-correlation of sur- 
face pressure with stopping muons rate, since a is large 
and negative. There is an anti-correlation of similar size 
with atmospheric height, as given by j3. Also, there is 
a weaker positive correlation with temperature at the 10 
kPa level, as given by 7 with a fairly large uncertainty. 

We tested the significance of the fit by repeating it with 
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FIG. 2. (Color online) Twenty-one weeks of stopping muon rate in 12 hour intervals expressed as a fractional change from the 
mean in the Model Fit Range. The error bars are statistical. The solid red curve is the result of a fit to the data in the Model 
Fit Range. 
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FIG. 3. (Color online) Detail of data from the Model Fit Range in Fig. [2] showing the rate dependence of the stopping muons 
together with the best-fit result to the model discussed in the text. The error bars are statistical. 



fewer parameters. The result with the temperature fac- 
tor 7 excluded is shown in Table |H| The quality of the fit, 
as estimated by the reduced % 2 , is almost the same, so we 
can conclude that the temperature of the atmosphere at 
the altitude of 10 kPa pressure is a statistically insignif- 
icant variable. Parameter a hardly changes, but /3 gets 
larger. There was scarcely any visual difference between 
the curves produced by these two fits. These resultant 
values are our best result for this measurement. 

To further test the model we restricted the parameters 



again to use only a for the ground-level pressure, exclud- 
ing the height parameter /3. This is shown in Table [TTT] 
The quality of the fit is now significantly worse, as seen in 
the reduced x 2 i an d the visual appearance of the fit was 
poor. Thus, ground level pressure alone is not enough to 
accurately track the rate of stopping muons. 

We used the same method to fit our data for the total 
event rate. This total rate at ground level in this detec- 
tor consists mostly of higher momentum muons (that do 
not stop), but with minor contributions from electrons, 
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FIG. 4. (Color online) For the same time period as in Fig. [5] the total cosmic particle rate expressed as a fractional change 
from the mean in the Model Fit Range. The error bars are purely statistical and generally smaller than the symbols. The solid 
blue curve is the result of a fit to the data within the Model Fit Range. The dashed red curve is the fit for the stopped muons 
from Fig. [3] 



TABLE III. Model parameters for the correlation of particle 
rates as in Table III but with only one free parameter. The 
third column is for muons stopping at ground level, while the 
fourth column is for the total detected cosmic particle rate. 



Parameter 


Eq 


tl 


Stopped 
Muons 


Total Particles 


Pressure 


a 




-3.3 ±0.5 


-2.0±0.1 








xl = 1-74 


X2= 2.22 



protons, and neutrons. The total rate includes the very 
small fraction (~ 0.4%) of stopped muon rate. The re- 
sultant fit is shown in Fig. |4j where one sees immediately 
the effect of the much higher statistics when counting all 
particles rather than just the lowest energy muons. The 
fractional change in total particle rate is now seen to vary 
quite smoothly on the timescale of days and weeks, al- 
beit with larger variations of order 10%. There are two 
model curves superimposed. The dashed red curve uses 
parameters from Table [I] obtained for the fit to stopped 
muons as per Fig. [3] The solid blue curve shows the new 
fit to the total particle rate with parameters also shown 
in Table [J] Within the Model Fit Range, the solid blue 
fit line is a close match to the data. Outside this range 
one sees again the large fluctuations we ascribe to space 
weather conditions that we could not track or reproduce. 
In fact, the fit to the regression model within the fit 
range led to a large value of xt of about 15. The statis- 
tics of this data set was so high that small fluctuations of 
unknown origin, presumed to be "space weather" , were 
evident. In order to compute legitimate uncertainties on 



the fit parameters in the least-squares algorithm, we in- 
creased the point-to-point uncertainties from 0.0010 to 
0.0035 to account for these fluctuations. This procedure 
did not affect the parameter values but allowed for the 
estimation of reasonable uncertainties for the correlation 
parameters. The fit result in the Model Fit Range is 
shown in Fig. [5j showing data points with their unmod- 
ified statistical uncertainties. The dashed blue curve fits 
the data very well. In this case the Pearson's correlation 
coefficient is 0.92, also indicating a very good fit. 

It can be seen that the amplitude of the fractional 
changes of the stopped muon rate (dashed red curve) 
is larger than the amplitude of the total particle rate 
changes (solid blue curve). However, they track each 
other very well over time. The fit parameters for the to- 
tal particle rate are included in Tables [TJ [ll]and III We 
see that a is about 50% larger for the stopped muons, 
while /3 is about four times as large. Taking into account 
the uncertainties on the values, a is larger for stopping 
muons at the "1.8 sigma" level, and j3 at the "5 sigma" 
level in Table [TTJ. From the figure, there can hardly be 
any doubt that the atmospheric dependencies indeed dif- 
fer between the total particle rate and the stopping muon 
rate. Thus we conclude that the lowest energy cosmic ray 
muons, the ones that stop in this detector, are several 
times more sensitive to ground-level pressure and atmo- 
spheric height variations than the total particle rate. 

It will be noted that our results do not include 
the ground level humidity measurements. Using local 
METAR data we found that the drastic variations in hu- 
midity from day to day and week to week were not cor- 
related in any detectable way to the stopping muon rate. 
Thus, we simply report that there was no measurable 
correlation. 
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FIG. 5. (Color online) Detail of data from the Model Fit Range in Fig. [i] showing the total cosmic particle rate expressed as 
a fractional change from the mean. The error bars are purely statistical. The curves are as in Fig. [4] 



Returning to results for the stopped unions, we as- 
sumed initially that the height and the temperature of 
the atmosphere at 10 kPa pressure was the optimal sin- 
gle level at which to determine the correlations. This was 
based on the knowledge [I] that this pressure level in the 
atmosphere is dominant in the formation of cosmic ray 
showers. However, one clearly expects the full develop- 
ment of muon showers to depend on the density of the 
atmosphere throughout an extended region, for which a 
full theoretical and numerical treatment is required [9]. 
Within the limited scope of our work we selected other 
reference altitudes from the IGRA data at which to com- 
pute the correlations for stopped muons to test the 10 
kPa assumption. The results are shown in Fig. [6] As 
can be seen, the xt values remain statistically close to 
unity except for the high and low extremes of altitude. 
This suggests that the parametrization we used works 
about equally well from ~ 5 to ~ 70 kPa, albeit with 
differing values for a, /3, and 7, wherein no single alti- 
tude seems more sensitive than the others. We believe 
that the fit worsens at the high altitude limit, below 10 
kPa, because the height and temperature there are in- 
creasingly above where the atmosphere has affected the 
primary cosmic particles. We believe that the increased 
uncertainties at the lowest altitudes arc due to a break- 
down of our assumption that the pressure measured in 
the laboratory is independent of the parameters mea- 
sured by the IGRA radiosondes. The model fails at low 
altitudes where the quantities are closely correlated. We 
note in the figure that the pressure parameter a tends 
to be smaller at lower altitudes. This is consistent with 
the trend reported in an early study for the total particle 
rate discussed in Ref. [IT] . 

In comparing the results of this measurement to previ- 
ous work, we have always converted the older results to 
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FIG. 6. (Color online) Evolution of the correlation parame- 
ters for stopped muons as a function of the reference pressure 
level (in kPa) used for atmospheric height and temperature. 
Solid red squares: ground-level pressure parameter a, solid 
blue triangles: atmospheric height j3, crosses: atmospheric 
temperature 7. The solid green diamonds are the reduced % 2 
values for the fits at each pressure level. 



our preferred unitless parameterizations for a, f3 and 7 as 
in Eq. [T] We used 16.6 km as the mean 10 kPa height of 
the atmosphere and 207 K as the mean temperature at 
this height. One of the previous experiments comparable 
to this one was performed by Trumpy and Trefall |12j . 
At ground level, they used stacks of counters separated 
by 10 or more centimeters of lead to record the total 
cosmic ray rates, and correlated these rates to the same 
atmospheric parameters as we do in these measurements. 
Table HVl includes their data for the total rate. One sees 
that their results agree with ours within errors for the 



TABLE IV. Comparison of the present results to selected previous measurements. All historic values were converted to 
consistent unitless coefficients. 



Parameter 


Eq.|l| 


Bernero 


Trumpy 


Chasson 


Shamos 


Cotton 


Grieder 


Fenton 






(This Paper) 


Ref.[E] 


Ref.pjB] 


Ref.[H] 


Ref.0 


Ref.pQ 


Ref.pl] 


Stopped Data 


Pressure 


a 


-3.2 ±0.5 


-3.5 ±0.3 


-3.6 ±0.5 


-7.1 








Altitude 


P 


-2.7 ±0.9 


0.0 ±0.3 


- 


- 








Temperature 


7 


+0.35 ±0.17 ±0.17 ±0.17 


- 


- 








Total Rate 


Pressure 


a 


-1.94 ±0.10 


-1.95 ±0.07 


-1.2 ±0.2 


-3.0 


-1.03 ±0.11 


-1.63 


-0.45 ±0.03 


Altitude 


P 


-0.8 ±0.2 


-0.53 ±0.08 


-0.54 ±0.08 


- 


-0.52 ±0.14 


-0.83 


-0.08 ±0.04 


Temperature 


7 


±0.08 ±0.04 ±0.10 ±0.04 ±0.14 ±0.04 


- 


±0.048 ± 0.056 


~0.2 


±0.04 ±0.02 



total particle rate. They then used a model to subtract 
the "hard" component, defined as tracks passing through 
at least 10 cm of lead, to arrive at a "soft" component. 
This soft component is not the same as our stopped muon 
measurement, but it is related. Additional "soft" contri- 
butions to their result could arise from electrons from 
in-flight muon decay, electrons from knock-on processes, 
and from very slow hadrons of all sorts. Nevertheless, 
we can compare their "soft" component with our results, 
and this is also shown in Table IIVI The values for the 
pressure correlation a are in agreement, but they were 
unable to obtain meaningful results for [3 and 7. 

Another of the scarce measurements of the "soft" cos- 
mic ray correlations was by Chasson [13j . Using stacks 
of Geiger counter telescopes with and without lead and 
iron absorbers between detectors, he factored a into 
"hard" and "soft" components. His result is included 
in Table [TV} For the total particle rate his results are in 
good agreement with ours, within uncertainties. For the 
stopped muon data he did not measure all three parame- 
ters, but gave only a value for the total barometer version 
of coefficient a. It is in agreement with ours. He did not 
measure stopping muons directly, as we have. 

We also mention the relevant results of Shamos and 
Liboff [14 . On the ocean surface they measured using a 
stack of ionization counters, so that, with some model as- 
sumptions, they factored a into "hard" and "soft" com- 
ponents. Their result is included in Table IV Unfor- 
tunately, they did not estimate uncertainties and they 
measured only the "total" coefficient, not the "partial" 
coefficients as we have. Thus, the poor agreement is dif- 
ficult to evaluate. Again, they did not measure stopping 
muons directly, as we have. 

Results from Cotton and Curtis [8] show poor to fair 
agreement with our results for the total particle rate. 
Their measurements were made at ground level using a 
stack of counters to select only the vertical component 
of the cosmic rate. Note that their temperature correla- 
tion is compatible with zero. Taking our experiment and 
these earlier experiments together, all give barely signif- 



icant positive values for 7; thus, there may be a hint 
for a small positive temperature correlation for the total 
particle rate at ground level. 

The text by Grieder [1] offers a set of total rate cor- 
relation parameters of unstated provenance and without 
uncertainty estimates which are in the same range as the 
other results in the table, including ours. 

For comparison of ground level and underground mea- 
surements we consider results from Fenton et al. |15j 
taken at 42 meters water equivalent depth, shown in the 
last column of Table |IV| For those muons, estimated to 
have at least 15 GeV energy, the correlations with mean 
production height and surface pressure are smaller than 
our results by roughly an order of magnitude, and their 
small temperature coefficient is consistent within errors. 
At greater depths it is has been shown [7] that the atmo- 
spheric effective temperature is the main correlator with 
high energy muon rate. Our small value for parameter 7 
at the surface is thus consistent with this picture. 



IV. CONCLUSIONS 

We have measured the rate of stopping muons at the 
surface of the earth and compared it with the rate of all 
charged particles. As was shown in Tables ||J [H] and |IH[ 
the atmospheric coefficients for the low energy stopping 
muons are significantly larger than for all particles. The 
linear regression model we have used represents both data 
sets well, in the sense that the fits to the data capture 
the behavior over a several-week period of low overall 
cosmic ray fluctuations, with good fidelity. The results 
are consistent with the picture that low energy muons, 
the least penetrating portion of the cosmic ray spectrum, 
are most sensitive to variations in atmospheric pressure, 
P, and the upper-level density related to H. We found 
no significant correlation with temperature at the mean 
production height, estimated using T. Our results were 
compared to previous measurements and found to be in 
fair to good agreement for the total particle rate. We 



were unable to find any previous measurements of the 
correlations for the stopping muon rate. 

It is also clear, as seen in Figures [2] and |4j that there are 
large and sudden departures from this model, which we 
ascribe to changes in the rate of incoming primary cosmic 
rays, albeit without independent evidence that this was 
the case. In the same figures it is seen that the model 
trends were sometimes qualitatively consistent with the 
rate changes in time periods not included in the fit region. 
However, certain time intervals exhibited unexplained 
step-like transitions in particle rates not accounted for 
in our phenomenological atmospheric weather data fit. 
Overall, we believe our most significant result is the set 



of correlation parameters for the rate of low energy muons 
that are stopped at ground level. We have shown that 
these correlations are significantly larger than those for 
the total rate of all cosmic ray particles. 
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